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' Abstract 

By using numerical simulation, we confirm that Takayasu-Sato-Takayasu (TST) model 
. which leads Pareto's law satisfies the detailed balance under Gibrat's law. In the simula- 

^ • tion, we take an exponential tent-shaped function as the growth rate distribution. We also 

^ , numerically confirm the reflection law equivalent to the equation which gives the Pareto in- 

Q ' dex /X in TST model. Moreover, we extend the model modifying the stochastic coefficient 

. under a Non-Gibrat's law. In this model, the detailed balance is also numerically observed. 

^ I The resultant pdf is power-law in the large scale Gibrat's law region, and is the log-normal 

distribution in the middle scale Non-Gibrat's one. These are accurately confirmed in the 



numerical simulation. 



C/3 

^ : 

PACS code : 89.65. Gh 



^ ! 1 Introduction 

O ■ Recent developments of computer technology enable vast quantities of economic data to be 

, analyzed. As one of remarkable examples, Fujiwara et al. [1] find that Pareto's law ^ 

^ '• P{x) oc x"(''+^) for x>xo (1) 

O 

can be derived from the law of detailed balance 

Pl2{xi,X2) = Pl2{x2,Xl) (2) 



and Gibrat's law [3] 

Q{R\xi) = Q{R) , (3) 

which are observed in massive amounts of economic data. Here x is wealth, income, profits, as- 
sets, sales, the number of employees and so forth, and xi, X2 are two successive those. Pi2{xi,X2) 
is a joint probability density function (pdf), and Q{R\xi) is a conditional pdf of the growth rate 
R = X2/X1 . 
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The Pareto's law ([T]) is not observed below some threshold xq O H], because the Gibrat's 
law dS]) does not hold in the middle scale region O [U [7] . In Ref . , we show that the detailed 
balance ([2]) and a Non- Gibrat's law 

QiR\xi) = Const. i?T*±(^i)-i ^ ^ 

t±(xi) = t±(xo)±a In— (5) 

lead the log-normal distribution in the middle scale region as follows: 

P{xi) = Cxi~(^+^) e'"^^ ^ for Xmin < Xi < Xq. (6) 

The Non-Gibrat's law means a statistical dependence of the growth rate distribution on the past 
value xi in the middle scale region. In the derivation, we employ the approximation (jj]) observed 
in profits data of Japanese firms [9], and find that the expression of t±{xi) ([5]) is unique under 
the detailed balance. Employing such empirical data, we confirm that the resultant distribution 
([6]) fits with data in the large scale Pareto's law region and those in the middle scale log-normal 
region consistently. The parameters are estimated as follows: a ~ for xi > xq, a ~ 0.14 for 
Xmin < xi < xolll Xo ~ 63, 000 thousand yen and 

Xmin ~ 

1, 600 thousand yen [8]. 
These data analyses are sufficiently rigorous, but at the same time are restricted in the finite 
period and category, because it is not easy to procure exhaustive data. If we build a model based 
on the detailed balance and (Non-)Gibrat's law, various economic situations can be simulated. 
This might make the reason clear why the parameters take the empirical values. Furthermore, 
we can study assets or sales data of firms which are difficult to obtain. In this paper, we will 
propose the simulation model based on the observation of economic data. 



2 TST model 

Firstly, we identify the model which leads Pareto's law based on the detailed balance and Gibrat's 
law in the large scale region. One of the simplest and powerful candidates is Takayasu-Sato- 
Takayasu (TST) model jtlO], which is given by the Langevin equation 

x{t + l) = b{t)x{t)+f{t), (7) 

where b{t) is a non-negative stochastic coefficient and f{t) is a random additive noise. They 
show that the conditions 

{\nb{t)) <0,{bitf)>l (8) 
are necessary and sufficient for the power-law ([1]), the index fi of which is given by the equation 

= 1 • (9) 

Here (• • • ) denotes an average over realizations. 

^ Here, a constant parameter a takes different values in two regions. This is not an exact procedure. However, 
for firms which are in the large scale region in both years {xi > xq and X2 > xq) or in the middle scale one 
(imin < xi < a:o aud Xmin < X2 < Xq) , this proccdurc is exact. In the database, most firms stay in the same 
region. This parameterization is, therefore, approximately valid for describing the probability density function 
P{xi). This is empirically confirmed in Ref. [8]. 
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In order to claim that this model is consistent with the empirical observation, the detailed 
balance and Gibrat's law must be satisfied. From Eq. ([7]), Gibrat's law is hold in the region 
where a noise f{t) is negligible. The result in Ref. [Ij suggests, therefore, that there is the 
detailed balance in TST model which leads Pareto's law. In the numerical simulation, we 
adopt the distribution of b{t) and f{t) to be Eqs. (j4])-((5|) with a = and Weibull distribu- 
tion: ^(^)'^~"'^ exp[— (^)'^], respectively. The typical scatter plot of the simulation is shown in 
Fig. [U where we take t+(xo) = 2.5, t-{xQ) = 1.5, xq = 10^, k = 0.5 and A = 30. To check the 
validity of the detailed balance, we take the one-dimensional Kolmogorov-Smirnov (K-S) test. 
We compare the distribution sample for P{xi G [lo3+0-2(n-i) ^ 10^+'^-^"], X2) with another sample 
for P{xi,X2 e [io3+0-2{"-i)^ io3+0-2n]) ^j^]^ n = 1, 2, • • • ,20 by making the null hypothesis that 
these two samples are taken from a same parent distribution. Each p value is shown in Fig. [2j 
The null hypothesis is not rejected in 5% significance level in the region where the noise f{t) is 
negligible. We recognize that the detailed balance is observed in the region. The lower bound 
of the region is denoted by Xmm, above which Pareto's law is observed (Fig. [3]). 
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Figure 1: The scatter plot of data points 
under Gibrat's law, the number of which is 
"500,000 ". 



Figure 2: Each p value of the one-dimensional 
Kolmogorov-Smirnov test for Fig. [TJ 




t+(xo) = a.o 

o-t+(X0) = 2.5 
o-i t+(XO) = 2.0 
M » tt(XO)-t-(XO) 



0.6 o.a 1 
T+(XO) -T -(XO) 



Figure 3: The cumulative distribution of 
P{x) shows the Pareto's law above Xmm ~ 
10^. The Pareto index /i is estimated by 
nearly 1. 



Figure 4: The relation between /x and 
t+(xo) — ^-(xo). Here, /i is estimated by the 
least square method of the simulation result, 
in which we take various values of t±(xo)- 



In addition, we also confirm the reflection law Q{R) = R ^ ^Q{R ^) derived in Ref. [Ij. In 
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this simulation, from Eqs.dH) and ([5]) with a = 0, the reflection law is reduced ton 

/i = t+(xo) - i-(xo), (10) 

which is equivalent to the condition Q. We observe this reflection law in various t±{xQ) (Fig.|3|). 
From these results, TST model is appropriate satisfying the detailed balance under the Gibrat's 
law. 



3 Simulation under Non-Gibrat's law 

If the detailed balance is satisfied even under a Non-Gibrat's law in the simulation, the log- 
normal distribution must be reduced in the middle scale region [8]. In a stochastic method, 
this scheme is not obvious analytically. Because a non-negative stochastic coefficient in Eq. ([7]) 
must be modified as b{x{t),t) under the Non-Gibrat's law. This extension exceeds the analytical 
framework of TST model. However, we can examine the scheme numerically. 

In the simulation, Xint is introduced satisfying t^{xint) = i~-(a^int)- We adopt the distribution 
of b to be Eqs. ([I])-® with a = for x > xq (Gibrat's law), a 7^ for Xjnt < x < xq (Non-Gibrat's 
law) and a = for x < Xint keeping t±{x) continuously (Fig. [5]). The last parameterization is 
imposed to exclude immoderate slopes of growth rate distributions. As the Non-Gibrat's region 
in the middle scale, we observe the region max(xmin, 3;int) < x < xq. 



t+(x) 




Figure 5: Continuous functions t±{x) = 
t±(xo) ± aln^ with a = for x > xq, 
a > for Xint < X < Xq and a = for 
X < Xint- The horizontal axis is logarithmic 
scale. 




Figure 6: The scatter plot of data points 
under Non-Gibrat's law, the number of 
which is "500,000 ". 



The typical scatter plot of the simulation is shown in Fig. [6l where we take t+(xo) = 2.5, 
t_(xo) = 1.5, Xq = 10^, k = 0.5, A = 30 and a = 0.1. To check the validity of the detailed 
balance, we take the same K-S test in the previous section. The result is shown in Fig. [71 The 
detailed balance is confirmed not only in the Gibrat's large scale region (x > xq = 10^) but also 
in the Non-Gibrat's middle scale one (max(xniini ^int 

) = 10"^ < X < xo = 10^). The resultant pdf 
of X is shown in Fig. [8l where the log-normal distribution in the middle scale region is observed 
in addition to Pareto's law in the large scale one. To confirm the validity of the log-normal 
distribution, we check the correlation between the parameter a inputted in the simulation and 



This equation is also confirmed empirically in profits data of Japanese firms [8]. 
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^ data point ^ 7 

Log-normal (10 <J< < 10"' 

— Pareto law (X > 10 ): p^O.976 




Figure 7: The result of one-dimensional 
Kolmogorov-Smirnov test for Fig. [6j 



Figure 8: The resultant pdf with t+(a;o) = 2.5, 
t-{xo) = 1.5, xq = 10^ k = 0.5, A = 30 and a = 0.1. 
In the large scale region, Pareto's law is observed. 
In the middle scale one, the log-normal distribution 
is observed. 



that estimated by fitting the pdf in the Non-Gibrat's region (denoted by agt). The result is 
shown in Fig. [9] for the case t+(xo) = 2.5, t-(xo) = 1.5, xq = 10^, k = 0.5 and A = 30 and in 
Fig. do] for the case t+(xo) = 3.5, t-{xo) = 2.0, xq = 10^, k = 0.5 and A = 30 for example. The 
correlation is very high. Consequently, this simulation model is also appropriate satisfying the 
detailed balance even under the Non-Gibrat's law. 



' ' ' 

0.1 0.2 0.3 

alpha 

Figure 9: The relation between a and agt 
for the case t+(xo) = 2.5, t-{xo) = 1.5, 
xo = 10^ /c = 0.5 and A = 30. 



^ ' ' ' 

0.1 0.2 0.3 

alpha 

Figure 10: The relation between a and a^t 
for the case t+{xo) = 3.5, t-(xo) = 2.0, 
Xq = 10^ k = 0.5 and A = 30. 



4 Summary and future problems 

In this study, by using numerical simulation, we confirm that TST model which leads Pareto's 
law satisfies the detailed balance under Gibrat's law. In the simulation, we take an exponential 
tent-shaped function as the growth rate distribution. We also numerically confirm the refiection 
law [Tj equivalent to the equation which gives the Pareto index fi in TST model. Moreover, we 
extend the model modifying the stochastic coefficient under a Non-Gibrat's law. In this model. 
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the detailed balance is also observed. The resultant pdf is power-law in the large scale Gibrat's 
law region, and is the log-normal distribution in the middle scale Non-Gibrat's one. These are 
accurately confirmed in the numerical simulation. 

In this simulation, we employ the Non-Gibrat's law ([5|) that the probability of the positive 
growth decreases and that of the negative growth increases symmetrically as the classification 
of X increases (Figs. [5l \TT\i . This is uniquely derived from the exponential tent-shaped form (j4]) 
observed in Japanese profits data [8j. On the other hand, as for sales or assets of firms, it is 
reported that the probability of the positive and negative growth decrease simultaneously as the 
classification of x increases (Fig. [T2|) . This difference is thought to be caused by the form of the 
exponential tent-shaped growth rate distribution. Eq. ^ is rewritten by using the pdf q{r\xi) 
for r = log^Q R as follows: 

^ogiQq{r\xi) = t±{xi) r for r^O. (11) 

For sales or assets of firms, these linear approximations are not used to express the growth rate 
distribution with the curvature [1]. 

A 



Figure 11: In the middle scale region, the ^. -.r, a r i r n 

ligure 12: As for sales or assets of firms, 

probability of the positive growth decreases , , , -,. r i • • t 

the probability of the positive and negative 

and the probability of the negative growth . 

growth decrease simultaneously as the clas- 

increases symmetrically as the classification 

sification of x increases. 

of x increases. 

In order to count the curvature, we add a second order term with respect to r: 

logiQq{r\xi) = c=ft±{xi)r + lnlOu±{xi)r^ for r^O. (12) 
In terms of Q{R\xi), these are expressed as 

QiR\xi) = d R-^Tt±{xi)+u±{x^)\nR R^l . (13) 

In this case, the detailed balance ([2]) is reduced to be 

-P(^l) ^ 1 Q{R~^\X2) ^ l+t+{xi)-t-{x2)-[u+{xi)-u-{x2)]lnR (-^^^ 

P{x2) R Q{R\xi) ^ ' 

for R> 1. By setting R = 1 after differentiating Eq. ()14p with respect to R, differential equations 
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are obtained. The solutions are uniquely fixed as 

t+ix) = 



i+(xo)+aln h — In — 

xo 2 



u^{x) 
u^{x) 
P{x) 



(q 

[3 + 5 



7]) In — 



Xo 3 
p-6 



^In^^ 



In2^-2ln3^ 
3 



Xo 



Xo 



X 



^+(2^0) - ^^T— In — - ^ In^ — 
Xo b Xo 



n_(xo) 
C x-^'^+i) 



- 2(5 , X 7 , 9 X 

In ^- In"^ — 

6 Xo 6 Xo 



exp 



-(a--) In — 
2 Xo 



2/? 



- In'^ In* — 

6 Xo 6 Xo 



(15) 
(16) 
(17) 
(18) 
(19) 



with r]/2 = — n+(xo) + n_(xo) . Same solutions are obtained for R < 1 . Here we set the 
parameters as follows: a,(3,^,6,r] = for x > xq (Gibrat's law), a,/3, 7, (5 / and 77 = for 
X < Xo (Non-Gibrat's law) because n±(xo) are constants. 

As an index of increase and decrease of the probability of the growth rate, we examine 
the differential coefficient at the origin: ^ logio '7('^|a^i)j,^_|_Q — +'^±{^1) ■ We are, therefore, 
interested in the increase and decrease of t±(x) in the middle scale region. In the case 7 > 
for example, the results are shown in Figs. [13] and [14] using the notation X = ln(x/xo), ^min = 
ln(xjnin/a;o) • Consequently the solution, that the probability of the positive and negative growth 
decrease simultaneously as the classification of x increases, exists in following parameter regions: 



a = and 5 < < and 7 = 0, 

a = and 5 — 7Xmin < (3 < and 7 > 

a = and 6 < (3 < —jX^ain and 7 < . 



(20) 
(21) 
(22) 



In any case, t±{x) have no first order term and the negative second order term with respect 
to ln(x/xo). Not only parameter /? but also 6 is necessary for the solution. We will apply the 
simulation to this analysis in the next work. 
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Figure 13: For 7 > 0, the monotonically Figure 14: For 7 > 0, the monotonically 

increasing or decreasing regions of t+(x) . increasing or decreasing regions of t_(x) . 

As other application of the simulation, the detailed balance is expected to be extended to 
the detailed quasi-balance [H] observed in the assessed value of land in Japan [121 [l3]- If 
these applications are studied adequately, we can manage the risk of these quantities by using 
knowledge of this method. 
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